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1 Introduction 

In 1955 Van Hove proposed a remarkable time rescaling procedure that 
enabled him to derive the master equation from the Schrodinger equation, for 
a quantum mechanical system endowed with an infinite number of degrees 
of freedom, such as a quantum field. The idea is to consider the limit 

A — * keeping i= \^t finite (A-independent constant), (1.1) 

where A is the coupling constant and t time. One then looks at the evolution 
of the quantum system as a function of the rescaled time t. This is called 
Van Hove's "A^t" limit and provided an interesting solution to some long- 
standing problems in quantum mechanics and quantum field theory, such as 
a rigorous justification of the Fermi "golden" rule and of the Weisskopf- 
Wigner approximation [^. 

Van Hove's prescription avoided the rigorous consequences of the quan- 
tum mechanical evolution law, which is governed by strictly unitary operators 
and predicts that the decay of an unstable quantum system cannot be purely 
exponential, being quadratic for very short times 0] and given by a power 
law for very long times [^ . These features of the quantum evolution are so 
well known that they are discussed even in textbooks of quantum mechanics 
P and quantum field theory [0 . The temporal behavior of quantum systems 
is reviewed in Ref. [§. 

One should notice that the problem of the deviations from exponential 
decay was considered an academic one until very recently. The renewed inter- 
est in the short-time nonexponential behavior was caused by a nice proposal 
by Cook Ip, the subsequent experiment performed by Itano et al |jlO| and 
the debate that followed [|lT|. One must notice, however, that Cook's idea 



and the subsequent papers did not deal with bona fide unstable systems. The 
latter require a quantum field theoretical analysis and the careful treatment 
of cut-offs and divergent quantities |T2|, ^ |T^. It is also worth emphasizing 
that no deviations from the exponential behavior for an unstable system were 
observed until 1997, when Raizen's group detected non-exponential leakage 



through a potential barrier |15 



In this paper we shall look at Van Hove's limit from the perspective of the 
complex energy plane. We shall consider some particular cases, concentrating 
our attention on the exponential decay law and the irreversible features [|16] of 
the evolution. There is interesting related work in the literature, in particular 

1 



in connection with quantum dynamical semigroups [O] and the so-called 



"stochastic limit" in quantum theory 18 



2 A simple example: A^-level atom 

We start our analysis by considering a simple nonrelativistic model: an A^- 



level atom in interaction with the electromagnetic field |jT9|]. This example 
will help us to pin down some salient features of the X^t limit. The Hamil- 
tonian is 

H = Ho + XV, (2.2) 



with {h = c = 1) 

Ho = Y. ^M^i^ + S / ^^ ^'^i 

/-oo p 



^ojl3^ojl3y 



Qlj/J 



(2.3) 
(2.4) 



M.i' /3 



where z/ runs over all the atomic states, 6j,, bu are the annihilation and creation 
operators of the atomic level z/, obeying anticommutation relations 



{bk.bl} = 6, 



U-, 



(2.5) 



and a^^, a^^|3 are the annihilation and creation operators of the electromag- 
netic field, satisfying commutation relations 



oP^ «!'/?'] =K^- ^0^/3/3' 



(2.6) 



where uj is energy and (3 stands for other (discrete) quantum numbers (e.g. 
13 = {j,m,e), where j is the total angular momentum (orbital+spin) of the 
photon, m its magnetic quantum number and e defines the photon parity 
P = (—1)-'"'"^'^^). The general features of the form factors are well known for 



a wide class of physical systems [^ and some particular cases of the above 
Hamiltonian have been widely investigated in the literature EI . 



Assume one can prepare, say at time t = 0, the system in the initial state 
\fi; 0) (atom in state /i and no photons). The problem of state preparation 
is a subtle one that will be carefully discussed later. The initial state is an 
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Figure 1: Graphic representation of (|2.10|) : S''^^ and S^^^ are in the first and 
second hne, respectively. 



eigenstate of the unperturbed Hamiltonian Hq and the evolution is governed 
by the unitary operator 



—iEt 

U{t) = expi-tHt) = f / dE^—, 

Ztt Jc h — H 



(2.7) 



where the path C is a straight horizontal line just above the real axis. By 
defining the resolvents {^E > 0) 



S{E)^{fx;0\ 



1 



|/i;0) 



1 



E -Un 



S\E)^{fi;0\ 



1 



E-H 



/i;0), 



(2.8) 



E — Ho 
Dyson's resummation reads 

s'{E) = s{E) + x^s{e)j:{e)s{e) + x^s{e)j:{e)s{e)j:{e)s{e) + ..., 

(2.9) 
where S(-E) = {fi;0\V{E — Hq) ^V\fi;0) is the 1-particle irreducible self- 
energy function, that can be evaluated by the expansion 



with 



S(E) = S(2)(E) + A'S(4)(E) + ... 






— UJ 



(2.10) 



(2.11) 



Both S*^^^ and S*^^) are shown as Feynman diagrams in Figure |l|. In the 
complex i?-plane S(£') and T,^'^\E) have a branch cut running from the 



ground-state energy to oo and no singularity on the first Rieniann sheet. 
Summing the series (|2.9|) one obtains 



We define the "survival" or nondecay amplitude and probability at time t 
(interaction picture) 

A{t) = (^;0|e*^«*f/(t)|/i;0), (2.13) 

Pit) = K^;0|e^^«*f/(t)|/i;0)p. (2.14) 

Incidentally, notice that the survival probability at short times behaves 
quadratically 

P{t) = 1 - t'/ri + ..., rz = (X'ifi; 0\V'\fi; 0))-'/'. (2.15) 

The quantity tz is the "Zeno time:" it is the convexity of P{t) in the origin. 
The nonexponential behavior at short times, besides its fundamental interest, 
entails the quantum Zeno effect [|]. Notice that the expansion ( p.l5|) is 



formal: we are implicitly requiring that the second moment of the interaction 



Hamiltonian exists — a delicate assumption in quantum field theory [|T2], |13 
The survival amplitude can be expressed as 

iEt 



iEtr"" •• I I ,„ e 



A(t) = — / dEe-'^'S'(E + uJo) = — / dE- -— — -. (2.16) 

In Van Hove's limit one looks at the evolution of the system over time inter- 
vals of order t = t/\^ {t independent of A), in the hmit of small A. Let us 
see how this procedure works in the complex-energy plane. To this end, by 
rescaling time t = A^t, we can write 

a(-L]=± [ dE-. ^-^ , (2.17) 

\Xy 27cJc E-E{X^E + ujo) 

where we are naturally led to introduce the rescaled energy E = E/X^. 
Taking the Van Hove limit we get 

T.{X^E + u^) ^ S(2)(A2e + o;^)^_ =S(2)(^^ + zO+) 

A — 



j^^p Jo UJ^-OJy-UJ + 20+ 2 



(2.18) 



where 

AW - V^Tdu ^£ML_, (2.19) 

r^) = 2'nY.W{uj)\\ (2.20) 

the term +iO+ being due to the fact that ^E > 0. The propagator becomes 
S'iE) = hm -^ J^ = — ^ (2.21) 

and the survival probabihty reads 

A{i) = \imA(^] =— I dEe-'^'S\E) = e-[iA{a.,)+r(a;^)/2]t^ (2.22) 
A^o \A ) 27r Jc 

which yields a purely exponential decay (Weisskopf-Wigner approximation 
and Fermi golden rule). In Figure ^we endeavoured to clarify the role played 
by the time-energy rescaling in the complex- E' plane. 

A few comments are in order. In the present model, the Van Hove limit 
works in two "steps." First, it constrains the evolution in a Tamm-Dancoff 



sector |2^: the system can only "explore" those states that are directly 
related to the initial state /x by the interaction V: the "excitation number" 
J\f^ = bjj^bf^+J2i3,oj 'Kjfs'^'^is becomes a conserved quantity and, as a consequence, 
the self-energy function consists only of a second order contribution that can 
be evaluated exactly. Second, it reduces this second order contribution, which 
depends on energy as in ( |2.11|) , to a constant (its value in the energy cj^ of the 



initial state), like in ( |2.18| ). Hence the analytical properties of the propagator, 
which had branch-cut singularities, reduce to those of a single complex pole, 
whose imaginary part (responsible for exponential decay) yields the Fermi 
golden rule, evaluated at second order of perturbation theory. 

Notice that it is the latter step (and not the former one) which is strictly 
necessary to obtain a dissipative behavior: Indeed, substitution of the pole 
value in the total self-energy function yields exponential decay, including, as 
is well known, higher-order corrections to the Fermi golden rule. On the other 
hand, the first step is very important when one is interested in computing 
the leading order corrections to the exponential behavior. To this purpose 




Figure 2: Singularities of the propagator 



T6p in the complex- i? plane. 
The logarithmic cut is due 

The 



The first Riemann sheet (I) is singularity free 

to T,^'^\E) and the pole is located on the second Riemann sheet (II) 
Van Hove rescaling procedure acts as a "magnifying glass" in the complex 
energy plane. After rescaling, the pole has coordinates ( p.l9| )-( p.20| ) in the 
complex-E" plane, without higher-order corrections in the coupling constant. 



one can solve the problem in a restricted Tamm-Duncoff sector of the total 
Hilbert space (i.e., in an eigenspace of A/"^ — in our case, Af^ = 1) and exactly 
evaluate the evolution of the system with its deviations from exponential law. 

3 A more general framework 

Let us generalize the analysis of the previous section. Consider the Hamilto- 
nian 

H = Ho + \V (3.1) 

and suppose that one can prepare an initial state \a) with the following 
properties 

iJok) = Ea\a), {a\V\a) = 0, 

{a\a) = 1. (3.2) 

The survival amplitude of state \a) reads 

A{t) = (a|e*^"*f/(t)|a) = ^ f dEe-'^'S'{E + Ea) 

= 2^ Ij^ E - \^i:{E + E,y ^^-^^ 

where S'{E) = {a\{E — H)^^\a) and S(i5) is the 1-particle irreducible self- 
energy function, that can be expressed by a perturbation expansion 

\^i:{E) = A2s(2) (E) + A^S(4) (E) + • ■ ■ . (3.4) 

The second order contribution has the general form 

1 ^^.. . ^.- 2 1 



S(2)(E) = {a\VP,- —P,V\a) = Y.\{a\V\n 

dE' V{E') 



E — Hq __l_ E — En 



2ti E-E'' 



(3.5) 



where P^ = 1 — \a){a\ is the projector over the decayed states, {|n)} is a 
complete set of eingenstates of Hq {HQ\n) = En\n) and we set Eq = 0) and 

T{E) = 27r ^ \{a\V\n)f6{E - E^). (3.6) 
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Notice that T{E) > for E > and is zero otherwise. In the Van Hove hmit 
we get 

where the resulting propagator in the rescaled energy E = E /\^ reads 

S'{E) = -, ^,J^ . .,,, , (3.8) 



E-S(2)(E, + 20+)' 

where we used 

i:{X^E + Ea) ^ S(2)(A2E + Ea 



A=0 



S(2)(E„ + ^0^ 



(3.9) 



(Weisskopf-Wigner approximation and Fermi golden rule). 

Let us compute the leading order corrections to the exponential behavior, 
in particular at short times. Just above the positive real axis we can write 



E(2)(^ + i0+) = A(E)--r(^), 



(3.10) 



where 



A(E) =v r 

Jo 



dE' T{E') 



2n E-E' 
We assume that T{E) is sommable in (0, +oo), so that for some r^ > 0, 

r(^) oc E''-^ for E^O. 



(3.11) 



It is then straightforward to obtain 



Tz 



te 
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/•oo 

Jo 
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-1/2 



A2r(E,) 



(3.12) 

(3.13) 
(3.14) 



which are the Zeno time and the lifetime, respectively. When time is rescaled 
according to Van Hove, the Zeno region vanishes 



rz = AVz = A 



'^ dE 

2^ 



T{E) 



-1/2 



0(A) 



(3.15) 



P(t) 



o(r) 




0(^^^"^) 



O(log(l/?i) ) 



Figure 3: Essential features (not in scale!) of the survival probability as a 
function of the rescaled time t. The Zeno time is 0(A), the lifetime 0(1), 
during the whole evolution there are oscillations of amplitude 0(A^''^^) and 
the transition to a power law occurs after a time 0(log(l/A)) [see ( p.l5| )- 
( 13.161) ]. The normalization factor becomes unity like 1 — O(A^). The dashed 
line is the exponential and the dotted line the power law. 



and the lifetime reads 



te 



AVe 



It goes without saying that the evolution must then be described in terms of 
the rescaled time i = \H. 

The details of the evolution were thoroughly investigated in |T^ in terms 
of the coupling constant. We only show in Figure ^the most salient features 
of the survival probability. In the Van Hove limit the coupling constant vanish 
(A -^ 0) and several things happen at once: the initial quadratic (quantum 
Zeno) region vanish, the oscillations are "squeezed" out and the power law 
is "pushed" to infinity: only a clean exponential law is left at all times, with 
the right normalization factors. All this is not surprising, being implied by 
the Weisskopf-Wigner approximation. However, the concomitance of these 
features is so remarkable that one cannot but wonder at the effectiveness of 



this limiting procedure. 

In atomic and molecular physics one gets very small deviations from the 
exponential law. For this reason, we displayed in Figure ^ the survival prob- 
ability by greatly exaggerating its most salient features. 

4 Relativistic quantum field theory 

We look now at a more complicated system. Consider the decay of a massive 
scalar particle $ of mass M into two identical massive scalar particles of 
mass m [See Fig. |5(a)]. We shall work in a completely relativistic framework. 
The Lagrangian density of our model reads 



\{d,^? -\mH^ + \{d,<Pf -\m^^ -^^ 



C = -{d,^f - -M2$2 + {d^<pf - -m202 _ ^^$02 ^ ^^^ (4_^^ 



The Lagrangian £ct contains the counterterms absorbing the infinite but un- 
observable shifts between the bare parameters (Mq, ttiq, Aq) and the physical 
ones (M, m, A): 

ccT = l^zid.^r - lhi<^' + lw,<pY - ls,n<p' - |/i$0^ (4.2) 

with 

Sz = Z -1, 6z = z -1, 
6m = M^Z-M\ 6m = mlz-m^, 6x = XozZ^/^ - X, (4.3) 

where Z and z are the field-strength renormalization constants ($o = 2'^/^$ 

and 00 = ^^/20) • 

The full two-point function 

G{p) = I d^x e^P-^(f]|T$(x)$(0) 1^) (4.4) 

is given by Dyson's resummation of the geometric series: 



p2 - M2 + iO+ p2 _ ]Vf2 + ^0+ ' '"^ '>2 _ M2 + iO+ 

i 
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(4.5) 



^ y 



p+k 



$ 



■z X ^^ 



\ 



' ^p 



(a) (b) 

Figure 4: (a) Decay of $ into (j) + (j). (b) Loop contribution to S*^^^(p^) 



where S(p^) is the 1-particle irreducible self-energy: 



(4.6) 



By using Feynman rules it is straightforward to write down the contribution 
of the loop diagram in Fig. Kb): 



(2) / 2\ _ izHJ^ f Cf^k 



^^looplP 



2 J (27r)4 P - m2 + zO+ (p + A;)2 - m2 + iO+ 



■ (4.7) 



^(2) 



We compute Sioop by using dimensional regularization. Introducing a Feyn- 



man parameter ^ and shifting the integration variable {q = k + xp) we get 

-'loop I 



-^^i±ip') 



2 7o 7 (27r)^[A;2-f2^fc-p + ^p2_ 



m 



212 



2 7o ^7 (27r)^[g2 + ^(l-^)p2_^2]2 



. (4.8) 



Performing a Wick rotation to Euclidean space (g^ = — ig°) and evaluating 
the momentum integral we obtain 



sSp(/) 



' r^ .. r d^qE 



-^1 d^ 



2 7o V (27r)^ [gi - ^(1 - ^p' + m^? 
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which diverges hke 2/(4 — D) for D ^ 4: 



^(2) 
''loop 



ip' 



/i^ 



2(47r)2 Jo 



di 



A-D 



7 — lo^ 






(4.10) 



where 7 is the Euler-Mascheroni constant. Let us define the particle mass 
M by the condition 

ReS(M2) = 0, (4.11) 

so that the propagator has the following behavior 



G{p) 



iZ 



p2_ M-^- zZImS(M2 



-, for p^ 



M^ 



(4.12) 



with Z ^ = 1 — S'(M2). Notice that Z is a finite field-strength normalization 
constant. Imposing the renormalization condition ( [4.11 ) one gets 



^^'\P') = S;Sp(P^)-ReSifip(M^ 



(2) 



/i- 



2(47r) 










where 



C 



d^log 



m^ -^{1-^)M^ 



m^ 



(4.13) 



(4.14) 



This is equivalent to setting 5m = ~ ^"^^^'^loipi-^'^) ^^'^ ^z = 0, in Eq. ( ^4.3|) . 
The function S*^2^(s) is analytic in the cut s plane and the discontinuity across 
the cut is (x > Am?) 



S(2)(x + zO+) - S(2)(x - zO+) = -27Tt-^Jl 



Am? 



X 



(4.15) 



Hence Tj^'^^s) can be represented by a dispersion relation. It is indeed 
straightforward to obtain 



S(2)(s) 



/i^ 



with 



2(47r)2 
P{s) 



Arrfi S'is — S'] 



Am? 



C 



(4.16) 



(4.17) 
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Therefore the propagator ( [4 .41 



has a simple pole Spoie near M^ in the second Riemann sheet and for s close 
to Spoie one gets 

G{s) ~ _'^ , (4.19) 

■5 'Spole 



where 



and 



Z = -f- = 1 + A2s(2)'(m2 + ^0+) + O(A^) (4.20) 

t ^nl'^poicj 



Spoie = M2 + Sn(spoie) = M2 + A2s(2)(M2 + z0+) + O(A^) 

= M'^-iX'^MT{M^) + 0{X^), (4.21) 

with 

,,2 

(4.22) 



The time evi 

At) -= 


Dlution of the o 
E G{t,p)=J 

J 2n 


orrelation function 

..-iEt 


(4.4) reads 
i 






E{2Ep + E)- 


- S(M2 + E{2Ep 


+ E)) 



(4.23) 

where Ep = y/p'^ + M'^ is the energy of the particle $. By introducing the 
rescaled time t = i/X^ and energy E = X^E, Eq. ([4.23|) becomes 

, f i \ jEp- r dE -fpT i 



^Xy J 27r E{2Ep + X^E)-^J:{M^ + X'^E{2Ep + X'^E)) 

(4.24) 
and taking Van Hove's limit, we obtain 

Ait) = lim e+^^U (^) = / ^ e-^^^G{E) (4.25) 
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with the hmiting propagator 



G{E) = hm- ^ 



A-o E{2Ej, + X^E) - ^S(Af2 + X^E{2Ep + X^E)) 
i 



2EpE-S(2)(M2+iO^ 
1 i 



(4.26) 



Therefore the time evolution ( |4.25| ) becomes 

T.~x 1 ( MT~\ 1 ( i \ , , 

and the particle decays exponentially with a mean lifetime 

r„ = :^r-i = ^~' , (4.28) 



M v^ 



f 



2' 



which has the proper relativistic time dilatation factor. This result is re- 
markably simple, for the model considered. 

5 Conclusions and comments 

The time evolution obtained in Van Hove's limit is always purely exponential: 
the quantum dynamics is governed by a master equation and by dynamical 
semigroups. However, it is obvious that the very procedure of time rescaling 
hides, in some sense, the problem of state preparation. 

In all the decay processes considered in this paper, an initial pure state 
is considered at time"t = 0." What is the meaning of t = 0? This question is 
often dismissed, in particular in quantum field theory, where all "relevant" 
physical quantities are constructed from the S'-matrix. There are however 
interesting examples in which the issue of state preparation is discussed, both 
in the context of semigroups [23| and scattering processes |24|. It is difficult 



not to wonder at the concept of initial time, in particular when one considers 
fundamental processes like particle creation in quantum field theory. Think 
again of the relativistic model analyzed in the previous section, as well as of 
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other examples recently considered in the literature [^, [T3|. Any classical 
picture of a decay process is necessarily mind-boggling. The preparation of 
an initial wave function (or initial state of a quantum field) is an inherently 
quantum mechanical process, certainly not an easy one to conceive. 

This problem is difficult to tackle. Nico Van Kampen, after refereeing one 
of our papers, put forward the following interesting and thought-provocative 



comment |26[: "As to your suggestion for preparing an initial pure state, 
there is no objection to it, from the mathematical viewpoint. But we are do- 
ing physics. Your construction is of the same calibre as the construction in 
statistical mechanics of those time reversed states whose entropies increase. 
The answer there too was that they are permissible from the mathematical 
point of view, but that they are tremendously improbable. Remember also 
the work by Wheeler and Feynman, which argued that there can be no co- 
herence in incoming waves owing to the absorbtion property of the universe. 
My feeling is that there is no real difficulty or paradox, but only the task to 
formulate precisely what is intuitively clear." 

We agree. Although an initial pure state like those considered in this 
paper are mathematically easy to conceive, their physical construction is pro- 
hibitive. Does this mean that nonexponential decays in atomic or elemen- 
tary particle physics are extremely improbable to observe because "unstable" 
quantum systems are practically always created in some sort of mixed states, 
whose statistical features justify time coarse-graining procedures like Van 
Hove's? This is an interesting question, which goes to the very core of the 
notion of irreversibility. A physical answer is needed. 
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